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I. INTRODUCTION 



As it was emphasized by Isern et al. |2j , due to simple cooling process in white dwarfs 
and precise measurements of luminosity curve became possible to use white dwarfs as a 
laboratory for analyzing problems of elementary particles physics. In the above works, it 
has particularly been suggested to study possible existence of axions on the basis of white 
dwarfs luminosity function. Following such an idea, we analyze the influence of lepton 
number violation on the luminosity of strongly magnetized iron white dwarfs. As is well 
known, the result of existence of Majorana type neutrino would be the lepton number 
violating process of electron capture by a nucleus X(A,Z) 

e~ + X(A,Z) -> X(A,Z-2) + e+ , (1.1) 

which is an analogue of the netrinoless double beta-decay, that is at present intensively 
studied [3|. 

Detailed study of strongly magnetized cold electrongas and its application to the magnetized 
white dwarfs has recently been done in Refs. [J 0]. This theory stems from Landau 
quantization of the motion of electrons in magnetic field 0, Q and of its modification to 
the case of very strong magnetic field It turns out that in systems with small number 
of Landau levels, which is restricted by a suitable choice of the strength of the magnetic 
field and of the maximum of the Fermi energy Ep of the electron gas, the mass of strongly 
magnetized white dwarf can be in the range (2.3-2.6) M , where M « 2x 10 33 g is the solar 
mass. It means that the strong magnetic field can enhance the energy of the electron gas to 
such a level that its pressure can force the gravity to allow the white dwarf topossess more 
mass than it is restricted by the Chandrasekhar-Landau limit of 1.44 M (9I. Ho|. 

Here we apply the above mentioned theory to study the influence of the double charge ex- 
change reaction (11.11) on cooling of the strongly magnetized iron white dwarfs. The threshold 
for this reaction with the initial nucleus ggFe an d the final one, IfCr, is A=6.33 MeV. So for 
the strongly magnetized iron white dwarfs with Ep > A this reaction can take place. We 
shall consider eF=EF/m e c 2 =20, 46, 90 and choose the strength of the magnetic field so that 
the value of the related parameter 7 will allow us to restrict ourselves to one Landau level. 

In Section [III we present methods and necessary input, needed for the calculations of the 
energy production due to the reaction (11.11) . including the double charge exchange width 
per one elementary reaction for chosen values of eF and 7. Then in Section [TTTl we give our 
estimate of the energy production and in Section IIV| we discuss our results and present our 
conclusions. At last, we provide in Appendix IA1 invariant functions, entering the positron- 
electron annihilation probability. Our main conclusion is that for sufficiently strongly mag- 
netized iron white dwarfs, the energy released in reaction (ll.ip can effectively retard their 
cooling. 



II. METHODS AND INPUTS 

In this section, we first discuss the necessary ingredients of the theory of strongly magnetized 
white dwarfs and then we present the formalism for the calculations of the double charge 
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exchange reaction (11.11) . 



A. Theory of strongly magnetized white dwarfs 



Theory of the strongly magnetized white dwarfs is based on Landau quantization of the 
motion of free electrons in a magnetic field B directed along the z-axis p, 0|. It is in detail 
discussed in Refs. Is}. In the relativistic case, one solves the Dirac equation, obtaining 
for the electron energy 



E v = m P c 



1 + — +2^ 



m t <-/ m\(? 



1/2 

(2.1) 



Here m e (e) is the electron mass (charge), c is the light velocity, p z is the electron momentum 
along the z-axis, v = I + 1/2 + a labels the Landau levels with the principal number I, 
a = ±1/2, and h is Planck's constant. The ground level (^=0) is obtained for 1=0 and 
er=-l/2, and it has degeneracy 1. Other Landau levels posses degeneracy 2. 

The last term at the r.h.s. of Eq. ( 12. ip can be written as 

ehB hjjH B 

2v— = 2v^- 2 =2v-. (2.2) 

micr m P c z B r 



e 



Here the cyclotron frequency uh = -^~ c and a critical magnetic field strength B c = = 
4.414 x 10 13 G. It is seen from Eqs. (12.11) and (12. 2p that the electron becomes relativistic if 
Huh > rn e c 2 , or if B > B c . 

In contrast to the density of electron states in the absence of the strong magnetic field, given 
as ( 27r )3fi3 d 3 p, the presence of such magnetic field modifies the number of electron states for 

a given level v to ^M^. dp z . Then the sum over the electron states in the presence of the 
strong magnetic field is given by 



E 



^(2^^/^ = (2^5:^/ d fe) ■ (2 - 3) 



Here 7 = B/B c , X e = h/m e c is the electron Compton wavelength, and g u = (2 — <5 0)l/ ) reflects 
the degeneracy of the Landau levels. 

The equation between the Fermi energy e-p and the Fermi momentum x-p{y) = pF(h>)/m e c 
for the Landau level, specified by v is obtained directly from Eq. (12.11) 

e 2 F =4(z/) + (l + 2z/ 7 ). (2.4) 

Then one obtains the electron number density from Eq. (12.31) as 

2 7 



n. 



(2tt) 2 A3 o 



^2 9v%f(v) ■ (2.5) 
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TABLE I: The values of the Fermi energy eF=EF/m e c 2 , used in the present study. Further n e - is 
the electron number density, p e - is the corresponding electron density, p m is the matter density, 
calculated according to Eq. (|2.7p for the nuclei 26 Ee. The values of 7 are the smallest values, 
satisfying Eq. (|2.6p and the mean energy e e - is defined in Eq. (|2.9p . 





n e -/10 33 [1/cm 3 ] 


p e -/10 6 [g/cm 3 ] 


p m /10 10 [g/cm 3 ] 


27 


e e - [MeV] 


20 


3.52 


3.20 


1.26 


400 


5.1 


46 


42.8 


13.4 


15.2 


2116 


11.75 


90 


321 


293 


114 


8100 


23 



The values of v in the sum are restricted by the condition that in Eq. (12.51) xf(u) > 0, which 
provides that e\ — (1 + 2b"~f) > 0, from which the inequality follows 

iW < integer - ) • (2.6) 

Supposing that the white dwarfs are electrically neutral, one deduces for the system of the 
electron gas and one sort of nuclei the matter density 

p m = fi e -mun e ~ = "^-vtla , (2.7) 



where \x e - = A/Z is the molecular weight per electron [A(Z) is the mass (atomic) number of 
the nucleus] , m\j is the atomic mass unit and uia is the mass of the nucleus with the mass 
number A. For the lightest nuclei, \i e - = 2, but, e.g., for j^Fe one obtains /i e - = 2.15. 

We also present the electron energy density at zero temperature 



2 7 m e c 2 — 



2 Vma 



7m e c 





{2-nYXl 

The energy per electron is then 



"max px p 

9u / [l+x 2 (z/)+2z/ 7 ] 1/2 dx(u 
Jo 

fmax 

Y.9 



x F (4 + 2,7) 1/2 + d + 2,7)/- F + (4 + 2Z/7)1/2 



(2.8) 



(l + 2z/ 7 ) 1/2 
e e -/n e -. (2.9) 



Numerically one gets e e - = Ep/2. 



Here we restrict ourselves to the strongly magnetized iron white dwarfs possessing the ground 
Landau level [y max = 0) and consider e F = 20, 46 and 90. Our choice of the parameter 7 
is 7=ep/2, which is the minimal value of 7, satisfying Eq. ( 12. 6ft . In Table HI we present, 
together with the values of e?, the related quantities, needed for the calculations. 
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B. Reaction rate 



The process of (e , e + ) conversion (II. ip is very similar to the neutrinoless double beta decay 
(Oz//3/3-decay) 

X(A, Z) -> X(A, Z + 2) + e" + e~ , (2.10) 

Both processes violate total lepton number by two units and therefore, take place if and only 
if neutrinos are Majorana particles with the non-zero mass. In addition, as we shall show 
below, the (e~, e + ) conversion rate is, like the Oi//3/3-decay rate, proportional to the squared 
absolute value of the effective mass of Majorana neutrinos \ (m u )\ 2 . This quantity takes the 
form 



Y^ulm- (2-11) 



i=i 



Here U is the 3x3 Pontecorvo-Maki-Nakagawa-Sakata unitary mixing matrix and m; (i = 
1,2,3) is the mass of the i-th light neutrino. Let us note that m v depends on neutrino 
oscillation parameters 6*i 2 , #i 3 , Am;? UN , Am^ TM , the lightest neutrino mass and the type of 
the neutrino mass spectrum (normal or inverted). 

From the most precise experiments on the search for the Oz//3/3-decay (ll-13j. by use of 



nuclear matrix elements (NME) of Ref. [14] . the following stringent bounds were inferred 

a 

\(m u )\ < (0.20 - 0.32) eV ( 76 Ge), 

< (0.33 - 0.46) eV ( 130 Te), 

< (0.17-0.30) eV ( 136 Xe). (2.12) 
However, there exists a claim of the observation of the 0^/9/3-decay of 76 Ge, made by some 



participants of the Heidelberg-Moscow collaboration [151 ] . Their estimated value of the 
effective Majorana mass (assuming a specific value for the NME) is |(?7v)| ~ 0.4 eV. In fu- 
ture experiments, CUORE ( 130 Te), EXO, KamLAND-Zen ( 136 Xe), MAJORANA/GERDA 
( 76 Ge), SuperNEMO ( 82 Se), SNO+ ( 150 Nd), and others O, a sensitivity 

\{m v )\ ~ a few 10~ 2 eV (2.13) 

is planned to be reached, which is the region of the inverted hierarchy of neutrino masses. 
In the case of the normal mass hierarchy, |(j7i„)| is too small in order to be probed in the 
0i//9/9-decay experiments of the next generation. 

For the sake of simplicity, the (e~, e + ) conversion on nuclei is considered only for the ground 
state to ground state transition, which is assumed to give the dominant contribution to the 
(e~,e + ) conversion rate. The spin and parity of initial ( 56 Fe) and final ( 56 Cr) nuclei in 
the ground state are equal, namely + . The incoming electron and outgoing positron are 
considered to be presumably in the Si/ 2 wave-states 



V>e- V2 Ve-) « y/F (Z,E e -)u(P e -), (2.14) 



^ e S + 1/2 \P e+ ) « y/F (Z-2,E e+ )u(P e+ ). (2.15) 
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The Coulomb interaction of electron and positron with the nucleus is taken by the relativistic 
Fermi functions F(Z,E e ~) and FiZ — 2,E e +) [l7j], respectively. Normalization of spinors 
is u^P)u(P) = 1 and v\P)v(P) = 1. P e - = (E e -,p e -) and P e + = (E e +,p e +) are the 
4-momenta of electron and positron, respectively. The above approximation is expected to 
work well for energy of incoming electron below 30, or even 50 MeV. 

The leading order (e~,e + ) conversion matrix element reads 

(f\S^\i) = 2rr5(E e+ - E,- + E f - Ei)(f\T^\i) , (2.16) 

with 

(f\T^) = i (m u )* ±-G}^F (Z, E e -WF (Z - 2, E e+ ) v(P e+ )(l + is)u{P e -) x 

4 M (^ + ). (2.17) 
R 

Here Gp = Gfcos9 c and E$ (Ef) is the energy of the initial (final) nuclear ground state. 
Later on, we neglect the kinetic energy of the final nucleus. The conventional normalization 
factor of NME, M^ el3+ \ involves the nuclear radius R — 1.2 A 1 ^ fm. For the weak axial 
coupling constant qa, we adopt the value qa = 1.269. 



The NME in Eq. fl2T6|) . defined as 



r (e/3+) 



M W+) = + M (^ + ) (2.18) 

contains the Fermi, M^ + \ and Gamow- Teller, Mq^. \ contributions. They take the fol- 
lowing form 

(e/3+) _ 4vr# r dq f 2 ( 2) 
Mp ~ (2tt)3 J 2q Mq)X 

4 r+e-^ln) (n\ £ m r+e**- r ™e^ r - |0.+) 

q-E e -+E n - Ei + ie n 

(0;i E m r+^\n)(n\ E/ T^+e-^-'-'e-^-'lO,^) ' 

q + E e+ +E n -E i + \e % ! ' ' 




= 4nR rdq 

GT (2tt) 3 J 2q JAKH ' 



(0/1 Ei T+g,e -"-'| n ) ■ (n| E m r+<r m e"'.- --e"'-"|0+) 
g — -E e - + — Ei + ie„ 



q + E e+ + E n - Ei + ie n * 

We use the conventional dipole parametrization for the nucleon form factors, normalized to 
unity 

M « 2) =( i+ £T 2 - /a( « 2 )=( i+ S"' (2 ' 21) 
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with My = 0.71 GeV, M A = 1.091 GeV. In the denominators of Eqs. (12191 and f[2~20D . E n 
and e n are the energy and width of the n-th intermediate nuclear state, respectively. 

For the considered (e~, e + ) conversion on 56 Fe, the typical momentum of intermediate neu- 
trinos is about 200 MeV (like in the case of the O^/3/3-decay [l8|), i.e., significantly larger 
in comparison with typical excitation energies of the intermediate nuclear states. Thus, in 
Eqs. (I2.19P and (12.201) . we complete the sum over the virtual intermediate nuclear states by 
closure, after replacing E n — Ei and e n with some average values {E n — Ei) and e, respectively 



|n)(n| 1 



q - E e - + E n - Ei + ie n q- E e - + (E n - E { ) + is ' 

\n)(n\ 1 
q + E e+ + E n - Ei + ie n ~ q + E e+ + (E n - Ei) + ie' 



As a result, the nuclear matrix element M„ e e \ decomposed into the contributions coming 
from direct and cross Feynman diagrams, takes the form 

= il4f +) + Mif a + \ (2.23) 

with 

M<f » - <0+, £ [ ,_ff^g* fa («< W, (-) 

M (e/3+) /n+ |V + +^ P io(qrim)f 2 (q 2 )qdg f 1 \ + 

M cVo - (0 >_. 7- r m — / — 07 • <x m - — 0J). 2.25 

j£ 71 Jo q + E e+ + (E n -Ei) +ie V Pi/ 

It is important to note that the value of = — E e -+ < E n — Ei > is negative for considered 
values of £7 e - and the studied nuclear system A = 56. Therefore, the contribution of direct 
Feynman diagram with the light intermediate neutrino has the pole at q = —E r — ie, as 
it follows from Eq. ( I2.24p . As a consequence, there is a non-zero imaginary part of the 
(e— ,e + ) conversion amplitude, which should be considered as well. 

The following comment is in order. In the expression (I2.19P for nuclear matrix elements 
we neglected the contributions of the higher order terms of the nucleon current 
(weak-magnetism, induced pseudoscalar coupling). As suggested by the analogy with 0uj3j3- 
decay, these terms should not be essential for the light neutrino exchange mechanism. 

Now we are ready to write down the expression for g.s. — > g.s. (e~ , e + ) conversion rate. The 
differential capture rate can be written as 

dT (e?+) = £ | < f\ T \i > | 2 2n5(E e - +E l -E } - E e+ ) , (2.26) 

where V is a volume of a phase factor. After the summed and squared T- matrix is calculated 

we get 



T^ +) = \(m u )\ 2 ^^(^=] F (Z,E e -)F (Z-2,E e+ )^ ' p e+ E e+ . (2.27) 



2 
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Here, E e + = E e - - A. 



The above result is for a single electron in a volume V. By assuming the density of electrons 
n e - (we replace 1/V with n e -), the reaction rate per nucleus is 



r ( ej 8+) 



m f 



m 2 167T 



1 (G p mj\ 
Stt 3 V y/2 J 



A 



(R 2 ml) 



F Q (Z,E e -) F (Z-2,E e 
2 p e + E e + n e - 



m 2 ml 



(2.28) 



For the reaction (11 .11) . occurring in the strongly magnetized white dwarfs, one should sum 
up in Eq. ( 12.281) over the electron energies according to Eq. ( 12.31) . This procedure provides 
the change 



F (Z, E e -) F (Z - 2, E e+ ) V -^- ^ _> 0( ep , 7 ) . 



The function 0(ep, 7) is defined as 

2 7 



(eF,7) 



m* m° e 



(e e - -Qf-l 



(2.29) 



(^fXlml J Q+l 
xF (Z,e e -) F (Z-2,e e+ )rfe e -, 



1/2 



(2.30) 



where e e ± = E e ±/m e c 2 and Q=A/m e c 2 . The calculated values of for chosen ep and 7 are 



0(20,200) = 1.80 x 10 3 , 
0(46, 1058) = 7.94 x 10 5 , 
0(90,4050) = 3.57 x 10 7 . 



(2.31) 



Then the reaction rate reads 



r (e/3+) 



m f 



Gpm 2 e x 



m 2 167T 3 \ y/2 ) [B?m\ 



M w+) 



(cf,7) • 



(2.32) 



Nuclear matrix element 



We use the Quasiparticle Random Phase Approximation (QRPA) [14j,ll8| to calculate nuclear 
matrix element for the transition (e~,e + ) on 56 Fe. For the A=56 system, the single-particle 
model space consisted of — 4froj oscillator shells, both for the protons and neutrons. The 
single particle energies are obtained by using a Coulomb-corrected Woods-Saxon potential. 
We derive the two-body G-matrix elements from the Charge Dependent Bonn one-boson 
exchange potential [19| within the Brueckner theory. The pairing interactions are adjusted 
to fit the empirical pairing gaps. 

The particle-particle and particle-hole channels of the G-matrix interaction of the nuclear 
Hamiltonian are renormalized by introducing the parameters g pp and g p h, respectively. The 
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calculations are carried out for g p h = 1.0. The particle-particle strength parameter g pp of 
the QRPA is fixed by the assumption that the matrix element M^ T of the lepton number 
conserving process of electron to positron conversion on nuclei with emission of neutrino and 
antineutrino 

e" + X(A, Z) ->■ X(A, Z - 2) + e + + u e + F e , (2.33) 

is within the range (0,0.30) MeV -1 . Recall that a comparable quantity Mq u t , associated 
with the two-neutrino double beta decay of 48 Se, 76 Ge, 82 Se, 128 > 130 Te and 136 Xe, does not 
exceed the above range by assuming the weak-axial coupling constant to be unquenched 
{(j a — 1-269) or quenched (gA = 1.0). 

As we already commented above, the matrix element ^ (see Eq. fl2.24p ) of the direct 

contribution contains an imaginary part that stems from the pole of the integrand at q — 
—E r — ie. The averaged energy of the intermediate nuclear states (E n — Ei) is assumed to 
be 5 MeV. Taking into account that the widths of low-lying nuclear states are negligible in 
comparison to their energies, one can separate the imaginary and real parts of this matrix 
element using the well-known formula 

= V--t7c5{a), (2.34) 

a + ie a 

valid in the limit e — y 0. 

In Fig. [H the absolute value of the nuclear matrix element M^ e/3+ ^ for 56 Fe is plotted as 
function of energy of incoming electron E e -. The allowed region is due to the uncertainty, 
associated with fixing the particle-particle parameter g pp . We found that the contribution 
from imaginary part of M* 6 ^ is small but increasing with E e - like the hole modulus of the 
(e~,e + ) nuclear matrix element. For further qualitative analysis of the (e~,e + ) capture rate 
we shall consider 

|M (e/3+) |^3, (2.35) 

for E e - from the range (6.33, 50) MeV. 



III. ENERGY PRODUCTION OF THE REACTION e -> e+ IN THE STRONGLY 
MAGNETIZED IRON WHITE DWARFS 



In order to estimate the energy production in one event, e r , in the reaction (11 TJ , we cal- 
culated the two-photon positron-electron annihilation probability per unique volume within 
the quantum electrodynamics framework [20|, 21] and integrated it over the energies of elec- 
trons €f = Ef/m e interacting with the positron in the final state of reaction fll.ip . according 
to the prescription (12. 3p . As a result, we obtained 

^ 2 7 



(2^ 



Jf 



0,1. 



In this section, we keep h 
Further, 

r (2t) 4 



c = 1 and our notations follow closely those of Ref. [20], ch. 8. 



(3.1) 



[ \M\ 2 5^( Pf +p e+ -k-k')f i (k + k' )dkdk', (3.2) 



€,€' . CF.CT' 
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FIG. 1: (Color online) The absolute value of nuclear matrix element M^ e/3+ ^ for 56 Fe as function 
of energy of incoming electron E e -. The allowed region of calculated nuclear matrix element is 
determined by fixing the particle-particle parameter to Mq^ from the range (0,0.30) MeV . 

fo{k° + k'°) = 1, fi(k° + k'°) = k° + k'°, the sum is performed over two photon linear 
polarizations e, e' and the average is done over the electron (positron) spin z-component a 
(cr'), the final photons have the 4-momenta k^ and k'^ (k° = \k\ = u, k'° = \k'\ = u/), the 
positron (electron) 4-momentum is p^ + (p^), also E e + = p° e+ and Ef = p®. In its turn, the 
amplitude M is 

e 2 

M = ~ A ,n ^ n — ^v( Pe+ ,(j'){4'[-i(ft f - Js)+m e ] 4/(Pf-k) 
4(2vr) 3 y/iuu 1 ) 

+ 4H(fo- fi')+rn e ] 4'/iPf-k')}uip f ,a). (3.3) 
The calculation of the \M\ 2 reduces to calculations of traces. As a result we obtain 
a 2 vr f +1 , u 



2m 2 e e e + J_! e f + e e + + (\p e +\ - \p f \)x 1 

+A 2 /((p f ■ k)(p f ■ k')) + A 3 /(p f ■ k'f } rniUu + u/) . (3.4) 

Here the bared quantities are expressed in the units of the electron mass and 

ml + E e+ E f -\p e+ \\p f \ ul = Ee++E (35) 
E e+ + E f + (\p e+ \ - \p f \)x ' } 

The scalar function A comes from the part of the traces that do not contain the factors 
(e • pf) and (e' • pf). The functions A iy i — 0, 1, 2, 3, are presented in the Appendix lAl 
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TABLE II: The values of the Fermi energy ep — Ep/ni e , used in the present study. Further R WD 
is the radius of the white dwarf, AL (AT) is the calculated change in the luminosity (surface 
temperature), and the released energy in single reaction e r is defined in Eq. (|3.6p . 





R WD /10 4 [m] 


AL [W/s] 


AT [K] 


e T [MeV] 


20 


42.3 


2.7 x 10 11 


27 


9.1 


46 


18.6 


3.3 x 10 14 


244 


25.2 


90 


9.4 


2.9 x 10 17 


1847 


49.1 



In the next step, we include Tj into the function 0(eF,7) , thus obtaining o (eF,7) and 
0i(ep,7) • The energy produced in one event per one second is then calculated as 

e r = 1 (e F ,7)/0o(e F ,7) • (3.6) 



In the interval of 1 year, the number of reactions in 1 cm 3 is n r = n m T, where n m is the 
number density of matter. Then the released energy in 1 cm 3 per 1 year is 

E = n m Te r ,[Jcm- 3 y- 1 ]. (3.7) 

Let us consider the white dwarf with the mass M WD =2 M Ri4 x 10 33 g. Its volume is 
V WD =M WD /p m , from which one obtains the radius R WD . In the full volume, the released 
energy per 1 year is EV [J y -1 ], from which one obtains directly the change in luminosity 
(released energy per 1 second), AL = EV/3.154 x 10 7 [W/s]. The influence on the surface 
temperature of the white dwarf can be calculated from the equation 



AL 



WD 



1/4 



Here s = 5.67 x 1CT 8 Wm^K" 4 . 

Next we present in Table [XT] the calculated change in the luminosity and in the surface 
temperature of the strongly magnetized iron white dwarfs, using the necessary input from 
Table H |(m„)|=0.3 eV and Eqs. (E5TJ , ( ES2D , (1235]) . Besides, R=1.2 A 1 / 3 « 4.59 fm. 

Our model of white dwarfs with the constant density is not fully realistic. On the other 
hand, more realistic model, considered in 0, IE[ for cf=20 shows (see Fig. 3b) that for the 
mass of the white dwarf M WD 2 M is R WD ~ 8.8 x 10 5 m, which is about 2 times larger 
than in our model. So in calculating the change in the temperature according to Eq. (13. 8ft . 
we used the radius twice as large as obtained from the model with the constant density. The 
correction to it for ep=46 and ep=90 can be estimated as follows. For ep=200, the value 
of the radius of the model 0, |H[ is R WD ~ 7x 10 4 m [22j, which is about 1.65 times larger 
than the constant density model R W d- Applying the linear interpolation to get a corrected 
value of R WD for 6f=46 (cf=90), we obtained that the calculated values of AT in Table ILT1 
are underestimated by ~ 1 % (4 %). 
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IV. DISCUSSION OF THE RESULTS AND CONCLUSIONS 



In this work, we studied the influence of the double charge exchange reaction fll.ip on 
the cooling of the strongly magnetized iron white dwarfs. This reaction is closely related 
to the neutrinoless double beta decay process, which is nowadays intensively studied (3)]. 
Both processes violate lepton number by two units and, therefore, take place if and only if 
neutrinos are Majorana particles with non-zero mass. For the case of light neutrino exchange, 
the (e~,e + ) conversion and the O^/3/3-decay rates are proportional to the squared absolute 
value of the effective mass of Majorana neutrinos, \(m u )\ 2 . 

As it is seen from our Table HI1 the double charge exchange reaction fll.ip can for sufficiently 
magnetized iron white dwarfs effectively retard their cooling by pumping the energy of the 
Fermi sea of electrons to the thermal energy of ions. It also follows that the above studied 
type of white dwarfs could, in principle, serve as a tool for observation of existence of the 
Majorana neutrino. The analysis of data on magnetic white dwarfs with B ~ 1-2 x 10 7 



G have very recently been published in Ref. |23|. As it can be seen from Table 4, error in 



provisional estimates of the temperature is at the level of 10 %. However, our estimates 
in change of the temperature due to the reaction (11.11) could be of interest for analysis of 
magnetized iron white dwarfs with B > B c = 4.414 x 10 13 G. 
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Appendix A: Calculations of traces 



Here we provide the invariant functions entering the positron-electron annihilation prob- 
ability (|3.4p . resulting from calculations of traces and summing over the photon linear po- 
larizations. The calculations are made in the Coulomb gauge, putting e =e'°=0 and using 

e i e j = ~ kikj , e i 6j = 5ij — fcjfcj . 



m A e A 



(k ■ k 



f\2 



- a 



mtA, 



mtAo 



(Al) 
(A2) 



m\A 3 
Further, 



(Pf ■ k )(Pf ■ k ') 

{(-a 3 + a 5 + a n )(p f -k) + [(p r k') - (k ■ k') + a t - a 2 )a 13 }/(p- kf 

= {[a x -a A + a 12 - a 13 + (a s - a 6 )/2](p f ■ k) + [a 3 - a w - a u + a 12 

+ (a 6 - a 8 )/2] (pj ■ k') + [— a\ - a 3 + a 4 + a n - a 12 + (a 5 + a 7 )/2)(k ■ k') 

+ (a 9 - a 3 )a u + (a 13 - a 4 )ai + a 3 a w } /(p s ■ k)(p f ■ k') , (A3) 

{[ (p f ■ k) - (k ■ k') + a 3 - a 9 ]a 10 + (- fll + a 4 + a 7 )(p f ■ k')}/(p f ■ k'f . (A4) 



a = 1 + (k ■ k'f , 

a l = (Pf ■ Pe+) ~ (Pf ■ k')(Pe+ ■ k') , 

«2 = (Pe+ ■ k) - (p e+ ■ k')(k ■ k') , 

«3 = (Pf ■ Pe+) ~ (Pf ■ k)(p e+ ■ k) , 

a 4 = (pf ■ k') - (p f ■ k)(k ■ k') , 
as = (Pf ■ Pe+) - (Pf ■ k)(p e+ ■ k) - (p f ■ k')(p e+ ■ k') + (p f ■ k)(p e+ ■ k')(k ■ kf) 

06 = (k- k') [-(py • jfe') + ( Pf -k)(k- k') 

a? = (Pf ■ Pe+) - (Pf ■ k)(p e + ■ k') - (p f ■ k')(p e+ ■ k') + (p f ■ k')(p e+ -k)(k- k') 



«8 



(k ■ k') 



a 9 



iPf ■ k) + ( Pf ■ k')(k ■ k') 
( Pe+ -k')-( Pe+ -k)(k-k'), 
aio = Pf - (Pf ■ k 'f , 
a n = [Pf ■ k) - (pf ■ k')(k ■ k') , 
ai2 =p}~ (P f ■ k) 2 - ( Pf ■ k'f + ( Pf ■ k)( Pf ■ k')(k ■ W) 

ais = Pf~ (Pf ■ kf ■ 

Here a = a/\a\ is the unit vector. The invariant function A arises from the part of traces 
that do not contain the factors (e ■ pj) and (e' • pj). If one puts Pf = 0, one obtains the 



(A5) 
(A6) 
(A7) 
(A8) 
(A9) 
(A10) 
(All) 
(A12) 
(A13) 

(A14) 
(A15) 
(A16) 
(A17) 
(A18) 
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positron-electron annihilation probability in the laboratory frame of reference. At threshold 
positron energies one gets for T 



a 2 ix 
m 2 e 



(A19) 



which provides for the annihilation cross section 

\ • (A20) 



a 2 ir 



vm c e 



where v is the positron velocity. 
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